Slime mould solves maze in one pass . . . assisted 
by gradient of chemo-attractants 
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Abstract 

^n^J Plasmodium of Physarum polycephalum is a large cell, visible by un- 

aided eye, which exhibits sophisticated patterns of foraging behaviour. 
The Plasmodium's behaviour is well interpreted in terms of computation, 
where data are spatially extended configurations of nutrients and obsta- 

Ph cles, and results of computation are networks of protoplasmic tubes formed 

by the Plasmodium. In laboratory experiments and numerical simulation 
we show that if Plasmodium of P. polycephalum is inoculated in a maze's 
peripheral channel and an oat flake (source of attractants) in a the maze's 

1-1 central chamber then the Plasmodium grows toward target oat flake and 

connects the flake with the site of original inoculation with a pronounced 
protoplasmic tube. The protoplasmic tube represents a path in the maze. 
The Plasmodium solves maze in one pass because it is assisted by a gra- 

ly-^ dient of chemo-attractants propagating from the target oat flake. 
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1 Introduction 



A typical strategy for a maze-solving is to explore all possible passages, while 
marking visited parts, till the exit or a central chamber is found. This is what 
Shannon's electromagnetic mouse Theseus was doing in first ever laboratory 



^ experiment on solving maze by physical means 11 . The task of maze search is 

time consuming for a single mobile computing device. Therefore with advent of 
unconventional computing paradigm scientists focused on uncovering physical, 
chemical and biological substrates which can solve maze in parallel. Thus in [To] 
it is experimentally demonstrated that detecting superposition of excitation 
wave fronts propagating from source to destination, and from destination to 
source, in a maze filled with Belousov-Zhabotinsky (BZ) reaction mixture allows 
us to approximate a shortest path in the maze. Travelling excitation waves 
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propagate to all channels of the maze thus examining the maze's structure in 
full. The BZ maze solver does not represent path by the medium's physical 
characteristics — an external observer is required to reconstruct a shortest path 
from chemical waves' dynamics. 

In a gas-discharge maze-solver [9 an electric field, generated between source 
and destination electrodes, explores a maze in parallel. The field is stronger 
along the shortest path between source- and destination-electrodes therefore 
channels of the shortest path glow with high intensity. Such maze-solver has a 
volatile memory because the shortest path is visible as long as voltage to the 
electrodes is applied. A slime mould maze-solver reported in [8 is a fusion of 
BZ and gas-discharge approaches. Plasmodium of P. polycephalum is placed 
in several sites of maze at once. Initially the Plasmodium develops a network 
of protoplasmic tubes spanning all channels of the maze, thus representing all 
possible solutions. Then oat flakes are placed in source and destination sites 
and the Plasmodium enhances the tube connecting source and destination along 
the shortest path. A selection of shortest protoplasmic tube is implemented via 
interaction of propagating bio-chemical, electric potential and contractile waves 
in Plasmodium's body. Not shortest or cul-de-sac branches are abandoned. Thus 
we can call prototype 8 as 'pruning Plasmodium'. The largest protoplasmic 
tube (exactly its walls), which represents the shortest path, remains visible 
even when Plasmodium ceases functioning. Therefore we can consider the slime 
mould maze-solver as having non-volatile memory (analogous to precipitating 
reaction-diffusion chemical processors). 

In laboratory prototypes of a mobile droplet p] and hot ice computer n] a 
computation of shortest path in a maze is separated on two stages: computation 
of many-sources-one-destination set of paths and extraction of a shortest path 
from the set of paths. Many-sources-one-destination paths are computed by 
spreading acidity in 6 and propagating crystallisation pattern in [l] . A shortest 
path is selected and traced by a droplet travelling along pH gradient [6] or a 
virtual robot traversing crystallisation pattern 1 . In present paper we apply 
the two-stage computation in design of slime mould maze-solver which navigates 
a maze in one go, without exploring all possible options because the options are 
'explored' by chemo-attractants emitted by destination oat flake. 

2 Experimental maze-solving 

In laboratory experiments we used plastic mazes (Tesco's Toy Mazes, Tesco 
Pic), 70 mm diameter with 4 mm wide and 3 mm deep channels (Fig. |TJl) . 
We filled channels with 2% agar gel (Select agar, Sigma Aldrich) as a non- 
nutrient substrate. We smeared top of channel walls with strawberry flavoured 
Chasptick (Pfizer Consumer Healthcare Ltd) to deter Plasmodium from making 
'illegal' shortcuts over the walls separating channels. A rolled oat was placed in 
the central chamber of the maze and an oat flake colonised by Plasmodium of 
P. polycephalum was placed in the most peripheral channel of the maze. Mazes 
with plasmodium were kept in the dark in 20-23 C° temperature. Images of 
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Figure 1: Experimental maze-solving with Plasmodium of P. polycephalum : 
(a) maze used in experiments on slime mould maze-solving; (be) Plasmodium 
is inoculated in peripheral channel, east part of the maze, and a virgin oat 
flake is place in central chamber; (b) scanned image of the experimental maze, 
protoplasmic tubes are yellow (light gray); (c) binarised, based on red and green 
components, images, major protoplasmic tubes are thick black lines; (d) scheme 
of Plasmodium propagation, arrows symbolise velocity vectors of propagating 
active zone. See experimental laboratory videos at |http : //www . youtube . comTl 
user/PhysarumMachines 

mazes were scanned in Epson Perfection 4490. Photos are taken using FujiPix 
6000 camera. 
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A typical experiment is illustrated in Fig. [T] We placed an oat flake in the 
central chamber and inoculated Plasmodium of P. polycephalum in a peripheral 
channels. The Plasmodium started exploring its vicinity and at first generated 
two active zones propagating clock- and contra-clockwise. By the time diffus- 
ing chemo-attractants reached distant channels one of the active zone already 
became dominant and suppressed another active zone. In example shown in 
Fig. []Jib active zone travelling contra-clockwise dominated and 'extinguished' 
active zone propagating clockwise. The dominating active zone then followed 
gradient of chemo-attractants inside the maze, navigated along intersections of 
the maze's channels and solved the maze by entering its central chamber. 



3 Mechanics of maze-solving: Numerical mod- 
elling 

A plasmodium of P. polycephalum can be seen as a network of coupled bio- 



chemical oscillators 12, 14 . Interactions between the oscillators determine 
space-time dynamics of contractile activity and protoplasmic streaming in pro- 
toplasmic tubes [3j[5] and ultimately shape of Plasmodium's cell [7]. Laboratory 
experiments show that oscillators in plasmodium network interact similarly to 
neurons in a simple neural networks 13 . We can speculate that Plasmod- 
ium's active zones (analogs of growth cones of maturing neuroblast) establish 
mutually inhibiting relationships. An active zone proximal to a source of chemo- 
attractants sends biochemical and electrical signals which suppress activity of 
active zones distal to the source. 

A profile of Plasmodium's active zone on a non-nutrient substrate is isomor- 
phic to shapes of and behaves analogously to wave-fragments in sub-excitable 
media [5] . When active zone propagates two processes occur simultaneously — 
movement of the wave-shaped tip of the pseudopodium and formation of a trail 
of protoplasmic tubes. We simulate the tactic traveling of plasmodium growth 
front using two- variable Oregonator equation |4j: 

du x 1 , o / n I \ ^X Q \ r-v 2 

-^r = -i u x ~u x - [fv x + <p x ) — ) + D u \/ u 

at e u x + q 

dv x _ 

dt ~ Ux Vx ' 

The variable u x is abstracted as a local density of Plasmodium's protoplasm 
at site x and v x reflects local concentration of metabolites and nutrients. Pa- 
rameters q and / are inherited from model of Belousov-Zhabotinsky medium. 
Parameter <p x characterises excitability of medium in the Oregonator model and 
can be seen as analog of a degree of Plasmodium's isometric tension response 
to chemo-attractants' concentration [16]. We integrate the system using Euler 
method with five-node Laplacian, time step At = 5- 10 -3 and grid point spacing 
Ax = 0.25, e = 0.03, / = 1.4, q = 0.022. Walls of the maze are impassable 
obstacles, grid nodes occupied by the walls are non-excitable. 
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Figure 2: Simulating maze-solving in Oregonator model: (a) gradient of chemo- 
attractants originated in central chamber, intensity of black is proportional 
to concentration of the attractants; impassable walls of the maze are white; 
(b) time-lapse images of propagating active zones, active zones are recorded 
every 400th step of simulation as sites with u x > 0.1 are shown by red 
(dark gray), impassable walls of the maze are green (light gray); site of ini- 
tial medium's perturbation, analog of Plasmodium inoculation, is indicated by 
arrow; (c) major network of protoplasmic tubes extracted from (b). See videos 
at http : //www . youtube . com/use r/PhysarumMachines 



Let c x be a concentration of chemo-attractant at site x. At moment t = 
only site occupied by oat flake has concentration c x = 1, all other sites have zero 
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concentration. At every moment t of simulation a site x updates its state c x by 
the following rule. If c x = and there is at least one immediate neighbour y such 
that c y > then c x +1 =t~ 1 . Such rude approximation of diffusion is enough to 
build satisfactory gradient of chemo-attr act ants for guiding propagating wave- 
fronts. The gradient developed in simulation is shown in Fig. [2^i. 

Initially all sites have the same excitability parameter <j) x = 0.04515. The 
less is value of <j) x the higher is excitability of x. The excitability is revised as 
follows. At every crth step of simulation we detect excited, i.e. u z > 0.1, site 
z with maximum, amongst all sites excited at this moment, concentration of 
chemo-attractants c z . Then for every site x we revise its excitability as follows: 
if c x < c z then <j) x = 0.09. In simulation illustrated in Fig. a = 500. The 
medium is perturbed by an initial excitation, where a 5 x 5 sites are assigned 
u = 1.0 each (Fig.J^). The perturbation generates a propagating wave-fragment 
travelling along gradient of chemo-attractants. Some travelling waves branch 
but branches are extinguished due to increase of sensitivity threshold. 

To imitate formation of protoplasmic tubes we store values of u in matrix 
L, which is processed at the end of simulation. For any site x and time step t 
if u x > 0.1 and L x = then L x = 1. The matrix L represents time lapse su- 
perposition of propagating wave- fronts. The simulation is considered completed 
when propagating pattern reaches destination site (central chamber in scenario 
illustrated in Fig. [2| and halts any further motion. At the end of simulation 
we repeatedly apply the erosion operation [2] to L. This operation symbolises a 
stretch-activation effect 5 necessary for formation of Plasmodium tubes. The 
resultant protoplasmic network provides a good phenomenological match for 
networks recorded in laboratory experiments (Fig. [2b). 



4 Discussion 

We experimentally demonstrated that Plasmodium of P. polycephalum solves 
maze in one-pass, i.e. without exploring all possible solutions, if a source of 
chemo-attractants is placed at destination site. In our experiments we inocu- 
lated slime mould in a peripheral channel and a target oat flake in central cham- 
ber. Positions do not matter, if there is an obstacle-free pass from a source to a 
destination the slime mould will trace it. For example, slime mould inoculated 
in a central chamber finds exit out of a maze. In Fig. [3] we placed an oat flake 
colonised by Plasmodium in the central chamber and put an attracting oat flake 
in the west part of an outer channel. The Plasmodium navigated thought first 
junction. It branched at the second junction. Active zone travelling clock-wise 
detected higher concentration of chemo-attractants than concentration detected 
by contra- clockwise propagating active zone. Thus active zone moving clock- 
wise became dominating. Eventually it reached its destination site (oat flake) 
along the shortest path (Fig. [3]). 

Is our prototype better — in terms of complexity or real-life speed and 
costs — then other laboratory prototypes of maze- solvers (Tab. [if? All existing 
prototypes have the same computational time complexity O(L), where L is a 



6 



(a) 



(b) 



Figure 3: Experimental maze-solving with Plasmodium of P. polycephalum 
: Plasmodium is inoculated in central chamber and a virgin flake is place in the 
most peripheral channel, west of the maze, (a) photo of experimental maze, 
(b) binarised image. 



Table 1: Brief comparison of laboratory prototypes of maze-solvers implemented 
in spatially extended physical, chemical or biological media. Two stages of 
maze-solving are outlined: maze exploration, or computation of many-sources- 
one-destination paths, and path tracing, or following from source site to specified 
destination site. 

Prototype Maze is explored by Path is traced by 



BZ medium [10] 
Pruning Plasmodium [7] 
Gas- discharge [9] 

Hot ice rn 

Mobile droplet 6 
Single-pass Plasmodium 



Excitation waves 
Plasmodium 
Electrical field 
Crystallisation pattern 
Diffusing chemicals 
Diffusing chemo-attractants 



Computer 

Plasmodium 

Electrical field 

Computer 

Droplet 

Plasmodium 



length of a worst-case scenario shortest path in a maze (a worst case is when 
shortest path spans all channels of a maze, as e.g. Chartres and Remis mazes 
shown in Fig. |4|. 

The time complexity is determined by time taken by waves, electricity, Plas- 
modium or diffusing chemicals to explore the maze. Al these substances explore 
maze in parallel by propagating from their original locations to all other sites 
of maze, thus it takes them O(L) steps to span the maze. Tracing of a path 
from source to destination is rather a menial tasks — to follow gradients or 
other physical or chemical characteristics of the medium. It takes O(L) steps 
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Figure 4: Solving Chartres (a) and Reims (b) mazes in Oregonator model 
of slime mould. Time-lapse images of an active zones (analogs of excitation 
wave-fronts in sub-excitable chemical medium) following a gradient of chemo- 
attractants emitted by central chamber of mazes. At the beginning of experi- 
ments simulated Plasmodium is inoculated at the mazes' entrances, south edges 
of mazes. 

as well. Space complexity of all prototypes is the same O(L): it takes O(L) 
instances (molecules, charges, micro- volumes of mixture or protoplasm) to fill 
in all channels during exploration stage. 

In terms of a real time gas-discharge [9] is the fastest maze-solver, it takes 
just hundred of milliseconds to solve the maze. Hot ice computer p] is the 
second fastest one, a maze of 0.1 m in diameter can be solved in few seconds. 
Gas-discharge and hot-ice solvers are followed by BZ-medium [lO] and mobile- 
droplet p] solvers: solution time is measured in minutes and hours. Slime 
mould based solvers are the slowest ones: it usually takes them a couple of 
days to solve a maze like one we used in present experiments. The prototypes 
can be arranged in the following real-costs (which includes consumables and 
laboratory equipment) descending order: gas-discharge, mobile droplets, BZ- 
medium, hot ice, and slime mould. Gas-discharge maze solver is most expensive 
while P. polycephalum solver can run literally for free. 

Accuracy of slime mould computing is far from ideal. Situation illustrated 
in Fig. [I] is a typical one: slime computes almost shortest path, compare with 
the shortest path from source to destination computed in Oregonator model in 
Fig.[2j In some situations, e.g. the one shown in Fig. [5] Plasmodium chooses the 
longest path from source to destination. Said that in neither of 35 experiments 
we undertook plasmodium failed to solve the maze, it solves builds a path from 
source to destination. 
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Figure 5: Plasmodium of P. polycephalum chooses longest path to maze's 
central chamber: plasmodium is inoculated in peripheral channel, east part 
of the maze, and a virgin oat flake is place in central chamber; (a) scanned 
image of the experimental maze, protoplasmic tubes are yellow (light gray); 
(b) binarised, based on red and green components, images, major protoplasmic 
tubes are thick black lines. 

Our final 'disclaimer' is that neither of slime mould maze-solvers can suc- 
cessfully compete with existing conventional computer architectures. However, 
P. polycephalum is an ideal biological substrate which represents all essential 
features of reaction-diffusion chemical computers yet encapsulated in an elastic 
growing membrane 2 . It can be treated as a meso-scale prototype of future 
growing nano-scale circuits. We envisage slime mould maze- solvers can also con- 
tribute toward designs of parallel amorphous drug delivery systems, or smart 
needles, which grow towards the target tissue. 
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